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The equation of motion for the expectation value of a scalar quantum field does not have the local 
form that is commonly assumed in studies of inflationary cosmology. We have recently argued that 
the true, temporally non-local equation of motion does not possess a time-derivative expansion and 
that the conversion of inflaton energy into particles is not, in principle, described by the friction term 
estimated from linear response theory. Here, we use numerical methods to investigate whether this 
obstacle to deriving a local equation of motion is purely formal, or of some quantitative importance. 
Using a simple scalar-field model, we find that, although the non-equilibrium evolution can exhibit 
significant damping, this damping is not well described by the local equation of motion obtained from 
linear response theory. It is possible that linear response theory does not apply to the situation we 
study only because thermalization turns out to be slow, but we argue that that the large discrepancies 
■ we observe indicate a failure of the local approximation at a more fundamental level. 
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OS '. I. INTRODUCTION 

The dynamics of a classical inflaton field is commonly assumed to be governed by a local equation of motion of the 
| form 

m : .... 

+ 3iT0 + 77(0)0 + V e g (0)=O, (1.1) 

where the friction term rj(4>)4> represents the dissipative mechanism through which inflaton energy is converted into 
particles. In particular, in 'warm inflation' scenarios [ESS 01; this term is taken to be large enough for a significant 
fraction of the energy to be dissipated during inflation. Assuming that the classical inflaton is the expectation 
tIh ' value of quantum scalar field <!>, its equation of motion can be derived straightforwardly from a suitable quantum field 
Oh, theory, and it is not of the form Ql.ip . For example, if couples to a scalar quantum field x through an interaction 
i<j0 2 X 2 , one obtains a term g<fi(x 2 ) m the equation of motion, and the expectation value depends on through the 
effective mass of the x particles, to 2 (0) = to 2 + g0 2 . The field x here might be the quantum part of <& (that is, 
$ = + X, with — ($)) or it might be an independent field. In any case, the expectation value (x 2 ) is a highly 
, nonlocal quantity: it depends in a complicated way on the entire history of 0, in contrast to the local equation (|l.l|l , 
k^J ' where 0, and are all evaluated at a single time t. This situation is quite generic: a Yukawa coupling g'tpy^ip to a 
spin-i field tp produces a 0-dependent mass for the tp particles and a term g' (tptp) in the equation of motion, and in 
C$ ' general the equation of motion contains a sum of terms of this sort, but no friction term. 

The question naturally arises, whether the true equation of motion might be adequately approximated, under 
suitable circumstances, by a local equation similar to (|1.1|) . An apparently rather plausible argument suggests that 
this is indeed possible when changes sufficiently slowly with time, and supplies an estimate of the friction coefficient 
77(0). To evaluate (x 2 (t)} at a fixed time t, one expresses 4>{t') at some earlier time t' as 

0(f) = 0(i) + (t> - i)0(i) + l -{t' - tf4>{t) + .... (1.2) 

The evolution of x(t') is governed by a Hamiltonian H(t') = H(x, 4>{t')) that depends explicitly on time through 4>(t'). 
Inserting the expansion H1.2(l . we have 

H(t') = H( X ,<t>(t)) + AH(t'), (1.3) 

where the first term, say -ffo(x) = H(Xi4>(t)) is independent of t' . Assuming that the state of the x held is always 
close to thermal equilibrium, the time-dependent correction AH(t') can be treated using the standard methods of 
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linear response theory, which yields the approximation 

(X 2 (*)> » (X 2 )° q + i /* dt'([AH(t'), X 2 (t)]) eq (1-4) 

where the expectation values are taken in the equilibrium ensemble associated with iJo- Keeping just the first two 
terms in (|1.2|) . we find that AH(t') is proportional to 4>{t), and in this way we recover an equation of the form 
1)1.1(1 . Using this or equivalent methods, estimates of the friction coefficient were found in the 1980s for a single 
self-interacting scalar field 0, j an d more recently (mainly in the context of warm inflation) for more complicated 
models U Hfl El El E3 

Plausible though these calculations seem to be, they rely on an unproven assumption: namely, that the true non- 
equilibrium dynamics can legitimately be replaced by the linear response of a state of exact thermal equilibrium. It 
is therefore important to ask whether the same results can be recovered from the slow-evolution limit of the non- 
equilibrium dynamics. Recently we have addressed this question by exhibiting a sequence of approximations 
through which the exact non-equilibrium dynamics can be reduced to a set of local, coupled evolution equations for 
4>{t) and quantities characterizing the state of the particles into which <f> decays. To recover the local equation of 
motion 1(1.11) . it is necessary that the solution of these evolution equations should have an adiabatic expansion in 
powers of derivatives of (f>(t). Somewhat surprisingly, we found, as reviewed in section [H] below, that no such solution 
exists. Formally, this means that the equation of motion for cf>(t) is intrinsically non-local, and does not have a 
time-derivative expansion. In particular, the friction coefficient r)(<fi) does not exist. 

This formal result indicates that the local equation of motion is never strictly valid, but that does not necessarily 
mean that it cannot serve as a reasonably accurate approximation. Intuitively, one might still expect that (at least 
under favourable circumstances) dissipative effects, which certainly are a feature of the non-equilibrium evolution, 
should be fairly well represented by the friction term calculated in linear response theory. It is this issue that we 
attempt to address here by numerical methods. More specifically, one can envisage a further approximation to the 
non-equilibrium evolution equations in which self-energies are replaced with the values they would have in a state of 
thermal equilibrium. Then (see section ITT1 below ) the solution of the resulting equations does have a time-derivative 
expansion, and one recovers essentially the friction coefficient obtained from linear response theory. In what follows, 
we try to assess (within a particular model) whether the local equation of motion obtained in this way provides a 
good approximation to the non-equilibrium dynamics. 

Section [H] describes the non-equilibrium evolution equations and the adiabatic approximation to these equations 
that we solve numerically. In section ITTT1 we present solutions for the evolution from two different initial states: one 
is a state in which the adiabatic approximation predicts underdamped motion (section IIII All , while in the other 
it predicts overdamped motion ( section 1111 Bj l. In both cases, we find that, although although the non-equilibrium 
dynamics does yield damping, this damping is not well described by the adiabatic approximation. A brief account 
of the numerical methods used to obtain these results is given in section IIII CI Finally, section IIVI is a self-contained 
summary, during the course of which consider what might account for the failure of the adiabatic approximation in 
the situation we have studied. 



II. APPROXIMATE EVOLUTION EQUATIONS FOR A SCALAR-FIELD MODEL 

As in [3 El we work with a Minkowski-space theory. Extending our analysis to a Robertson- Walker spacetime 
presents no difficulties of principle, but it introduces complications which seem to have no bearing on the issue we 
address. Consider, then, a quantum field theory defined by the Lagrangian density 

C = - U($) + Idtf&X - k m U 2 - irAx 4 - £s*V- (2.1) 

It should be emphasized at once that we will not attempt to study the exact non-equilibrium dynamics of this model. 
In fact, we will arbitrarily discard various terms in its equations of motion which, though they may well affect the 
dynamics significantly, are irrelevant to the question on which we wish to focus. The first such approximation is to 
replace the quantum field Q(x,t) with its expectation value <j>(t), where we assume a spatially homogeneous state. 
Then the equation of motion for <fi(t) is 

+ U'(0)+^<X 2 OM)} = O (2.2) 
while the dynamics of the remaining quantum field x( x , t) is governed by the Lagrangian 



(2.3) 
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with m 2 (t) — m 2 , + g<j) 2 (t). An approximate treatment of this non-equilibrium dynamics yields, for the expectation 
value in 12. 2|) (which is independent of x for a homogeneous state) 

/d 3 k 
W)^) [1 + 2nk{t) + mei/kit)l {2A) 



where ui k = \/k 2 + m 2 (t) and the functions n k (t) and v k (t) : which characterize the state of the system of \ particles, 
are solutions of the evolution equations 

ed t n k (t) = a k (t)-T k (t)[l + 2n k {t)]+e^§-Reiy k (t) (2.5) 

ed t v k (t) = -2i[uj k (t)-iT k (t)]iy k (t)-a k (t)+e^^-[l + 2n k {t)], (2.6) 

In these equations, a k (t) and T k (t) are given by 

3 u S(u> kl + u k 2 - u) ks - uJ k )S(ki +k 2 - k 3 -k) 



a k {t) = J d 3 k 1 d 3 k 2 d 3 k 



(J k U) kl U) k2 L) k3 

x [(1 + n kl )(l + n k2 )n k3 + n kl n k2 (l + n ka )] (2.7) 

r fc ( t ) =°[ d^k^k^ ± ^ ^ UJk)s{kl ± fc2 ~ ** k) 

<±K l J U} k UJ kl UJ k2 UJ k3 

x [(1 + n kl )(l + n k2 )n k3 - n kl n k2 (l + n ka )] (2.8) 

with c = A 2 /64(27r) 3 , while e, which has the value e = 1, is a formal parameter introduced to facilitate an adiabatic 
expansion. 

The approximations that lead to these equations are described in detail in |15| . Essentially, they arise as a solution 
to the Dyson-Schwinger equations for 2-point functions, via a local Ansatz for self-energies, which are evaluated at 
2-loop order in perturbation theory. For a free-field theory, with A = 0, which also means a k (t) = T k (t) = 0, the 
evolution equations are exact (see, for example 0) and the functions n k (t) and v k {t) can be identified in terms of 
creation and annihilation operators for x particles: 

(4(t)a fc /(t)) = {2n) 3 S(k-k')n k (t) (2.9) 
(a k (t)a k ,(t)) = {2n) 3 5{k + k')is k {t). (2.10) 

With interactions included, l|2.5|l can be interpreted as a Boltzmann equation for number densities n k (t) ; the quantity 



S k (t) - a k (t) - T k (t)[l + 2n k (t)} (2.11) 
is the standard scattering integral associated with two-body elastic scattering, while the source 

m = ^Rev k (t) (2.12) 

represents particle creation and absorption due to the time-dependent mass m(t). 

Although the evolution equations (|2.5(l and Ij2.6|l are local, their solution is not. One can attempt to obtain a 
solution in the form of a time-derivative expansion, by writing 



M*) = n*(t) + e5n k (t) + 0(e 2 ) (2.13) 

eq 

v 



Mt) = 4 q (t)+eS,y k (t)+0(e 2 ), (2.14) 



expanding in powers of e and finally setting e = 1. However, it turns out that no such solution exists, for the following 
reason. To find the leading correction Sn k , one has to solve the integral equation 

/ dk'K(k, k')8n k , = Otn^ - ^Re (2.15) 
Jo 



where K(k, k') — {8S k /5n k /) e<t . It is straightforward to show that conservation of energy in two-body scattering 
implies the sum rule / °° dk k 2 Lo k K(k, fc') = 0. The right-hand side of (|2.15|l does not respect this sum rule, so the 
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equation is not self-consistent and the desired solution does not exist. (For elastic scattering, there is a second sum 
rule associated with particle- number conservation, which is also not respected by (I2.15|) : however, corrections to 
the Boltzmann equation (|2.5Jl from higher orders of perturbation theory would introduce inelastic scattering terms, 
destroying this second sum rule.) 

To the extent that our approximations afford an adequate description of non-equilibrium evolution (an important 
caveat, upon which we comment further in section llVjl this shows that the equation of motion (|2.2J) does not have a 
time-derivative expansion and cannot be represented by a local equation such as However, we can recover a time- 

derivative expansion by making a further approximation. The functions a k (t) and T k (t) are related to quasiparticle 
self-energies *£> k (t) (whose exact meanings are given in ^^|) by 

ak[t) = ~4^t) KW + ^Wl (2 - 16) 

r*(t) = -E-(t)]. (2.17) 

Let us replace these self-energies with constant equilibrium values. In what follows, we will call this the adiabatic 
approximation. Using this approximation, the time-derivative expansion can be carried through without difficulty. In 
particular, a k {t) and T k (t) no longer depend on 5n k (t), and the troublesome equation (|2.15(l does not arise. Since 
Co k {t) is proportional to <fi(t), so are 5n k (t) and 8v k {t), and these contributions to (x 2 ) give rise to a friction coefficient 



P(col + 3Tl)n k (l + n k ) , r k (u 2 k - 3rf)(l + 2n fc )" ,., |S( 



Here, n k stands for the equilibrium Bose-Einstein distribution n£ q = (eP Uh — l) , while T k is found by substituting 
this equilibrium distribution in (|2.8|l . This expression agrees, for example, with the result given in [7J in the limit 
T k << uj k considered by these authors. Roughly speaking, the second term represents a loss of inflaton energy arising 
from particle creation, while the first represents a transfer of energy to particles already present. 

This agreement is not surprising: the adiabatic approximation we have just introduced is roughly equivalent to 
the assumptions on which linear response theory is based. Intuitively, it would seem that this should be a fair 
approximation if (f>(t) does not change too rapidly, and if the state of the \ particles does not stray too far from 
an equilibrium distribution. We now attempt to test this intuition by comparing numerically the non-equilibrium 
evolution resulting from (|2.5|) and (|2.(i|) with that implied by the adiabatic approximation. 

To simplify matters a little, we observe that the frictional effect we want to study arises from the last two terms in 
the expression 1)2. 4f> for (x 2 ), and will discard the first term. Thus, we study the equation of motion 

4> + V'{<t>) + , 9 0<x 2 Oe, i))tru„c = (2.19) 

with the truncated expectation value 

(2^)32^) MO + Re > ( 2 - 20 ) 
the evolution of n k (t) and v k {t) being governed by l|2.5|l and 12.6fl . The adiabatic approximation yields the estimate 

g<f>(x 2 (x, t))£* = AU'(0) + n(4>)4 (2.21) 
where is given by (|2.18|l and the contribution from the leading order terms of the adiabatic expansion is 

2tt z j uj k 

with 

Re^ = - r 5* + 2 !f? . (2.23) 

In order to evaluate (x 2 ( x , 0)trunc> we nee d values for the inverse temperature f3 of the assumed equilibrium state. 
To this end, we observe that the evolution equations (|2.19|) . I|2.5|) and l|2.6|) conserve exactly the energy given by 

1 • 1 r°° 

E = -0 2 + V(0) + — J alkk 2 uj k n k . (2.24) 



5 



We will determine the evolving temperature in the adiabatic approximation by requiring this energy to remain 
constant. 

To avoid confusion, it may be worthwhile to reiterate that our goal is to test whether the expectation value (|2.20(l 
is adequately represented by the adiabatic expression l|2.21ll . We will do this by comparing solutions to the non- 
equilibrium evolution equations (|2.19() . (|2.20() . I|2.5[) and l|2.6|) with the adiabatic approximation to these equations, 
which consists of the equation of motion 

+ 77(0)0 + V'{4>) +AV\4) = 0, (2.25) 

with rj((f) and AV'(<p) given by H2.18fl and 12.2211 respectively, and (3 determined by energy conservation. The fact 
that we have arbitrarily discarded several terms does not affect this comparison, though of course it does mean that 
we are not in a good position to study the true dynamics of the original model l|2.1|l . Equally, our approximations 
are based in part on perturbation theory but the neglected higher-order terms do not affect the comparison, and we 
will not be concerned with choosing parameter values so as to make these terms small. 

III. NUMERICAL RESULTS 

We have obtained numerical solutions of the evolution equations presented in the previous section, taking the 
simplest available potential, namely 

V(0) = l -M 2 4> 2 . (3.1) 

This section presents two sets of solutions. In the first, the motion of is an underdamped oscillation; in the second, 
the motion is predicted by the adiabatic approximation to be overdamped. In each case, we solve, on the one hand 
the non-equilibrium equation of motion (|2.19|) together with l|2.5[) and (|2.6[) . and on the other hand the adiabatic 
approximation to this equation (|2.25|) . starting from identical initial conditions. At the end of the section, we comment 
briefly on the numerical methods used. 

A. Underdamped motion 

Figure n shows the motion of obtained in the two approximations, starting from a more or less arbitrary initial 
state. Throughout, the units are fixed by taking the mass of x particles to be m x = 1. Here, we chose the parameter 
values c = 0.5, g — 1.0 and M 2 = 0.4. The initial state for was = 5-0, = and the initial state of the \ 
particles was a state of thermal equilibrium with inverse temperature /3 = 0.3. With these values, the energy of 
the particles is about 86% of the total energy, so one might expect a significant transfer of further energy to these 
particles as well as a frictional effect from particle creation. The two superimposed oscillations in figure ^ are n °t 
easily distinguished, so we have plotted envelopes for the two curves: the solid envelope indicates the solution to the 
non-equilibrium evolution equations, while the dashed one corresponds to the adiabatic approximation. At first sight, 
there seems to be a similar rate of damping in each case, but there are significant differences. The friction coefficient 
of the adiabatic approximation given in l|2.18|l is a rather complicated function of 0, but suppose that it can be written 
as ti(0) = 2 ?}(0), where ?)(0) varies relatively slowly with 0. This suggests that large-amplitude oscillations should 
be strongly damped, while at small amplitudes damping would be ineffective. This behaviour is clearly apparent from 
figure^ However, the solution of the full evolution equations shows a different tendency: damping is smaller at large 
amplitudes than that predicted by the adiabatic approximation, and it does not disappear at small amplitudes. 

Thus, in this situation, although the non-equilibrium evolution equations do yield a damped motion, the damping 
is not well represented by the friction term as estimated from the adiabatic approximation. This is not necessarily 
unexpected, since both approximations yield fairly rapid oscillations, and the adiabatic approximation might not be 
expected to work well. We therefore turn to a different initial state which, according to the adiabatic approximation, 
ought to result in overdamped motion and slow evolution. 

B. Overdamped motion 

Suppose, as above, that the adiabatic friction coefficient can be represented as ?7(0) ~ t)4> 2 an( i al so that the 
adiabatic correction to the effective potential 12.22f) is times a slowly varying function. Then the adiabatic equation 
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FIG. 1: Underdamped motion of <j>(t) as calculated from the non-equilibrium evolution equations and via the adiabatic approx- 
imation. Envelopes of the two oscillatory curves are drawn to aid the eye. 



of motion (|2.25[l becomes 



+ + = 



(3.2) 



where t) and fj, 2 — M 2 + AV' /(f) are, at least over some period of time, approximately constant. If cf) can be neglected, 
then this equation has the overdamped solution 



m z - 



2 2 M 2 t 



1/2 



with a "terminal velocity" 



Under these circumstances, we have <f> w — /i 4 /?} 2 3 , so overdamped motion should occur if \(f>\ ^ \n 2 <j>\: or 

,2 



< 1. 



(3.3) 



(3.4) 



(3.5) 



Intuitively, the adiabatic approximation ought to be good if, in addition to this condition for neglecting 0, the state 
of the x particles thermalizes rapidly. According to the Boltzmann-like equation (|2.5|l . the mode of the particle 
distribution relaxes towards equilibrium with a characteristic relaxation time = 1/IV Taking the k = mode as 
representative, one can guess at a rough criterion for efficient thermalization, namely that the fractional change in 



the effective particle mass m 2 



+ gcj) 2 over a time Tfc =0 should be small. That is, Tk=o\dm /dt\ m 2 or 



< 1. 



fe=0 



(3.6) 



7 




FIG. 2: Motion of 4>{i) as calculated from the non-equilibrium evolution equations and via the adiabatic approximation. As 
determined from the adiabatic approximation, the initial conditions should lead to overdamped motion. 



Ideally, we would like to study a state in which both conditions (|3.5I) and (|3.(i[) are well satisfied, but this proves to 
be quite difficult. In particular, the second condition seems hard to achieve. Using i|3.4[l . it can be rewritten as 



< 1 



(3.7) 



and the thermal effects that make i) and Tk large also tend to produce large values of fi 2 . The dependence of ri((j)) 
(equation 1)2. 18fl ) on the parameters at our disposal is so complicated as to make a systematic survey of the parameter 
space impractical. By trial and error, we have chosen a set of values which yield 



5.9 x 1(T 



and 



2g4>4> 



m 2 r k=0 



0.97. 



(3.8) 



This means that according to the adiabatic approximation, the system is strongly overdamped, but thcrmalization 
is only moderately efficient. A careful, but by no means exhaustive, exploration suggests that this situation cannot 
be much improved within the model studied here. The parameter values leading to this state are c = g = 1.0, 
M 2 = 0.005, (3 = 0.09 and = 80.0. The terminal velocity O turns out to be = -0.08967 . . .. 

Taking these initial conditions, we find for <fi(t) the motion depicted in figure Evidently, the motion obtained 
from the adiabatic approximation (the thinner curve) is indeed overdamped at early times. In fact the solution agrees 
well with (|3.3() for times up to about t = 150. But as <f> becomes smaller, the condition (|3.5() for overdamping no 
longer holds and, as in figure ^ the motion at small amplitudes is essentially undamped. 

The solution of the non-equilibrium evolution equations (the thicker curve in figure^ is strikingly different: although 
there is significant damping, the motion is far from being overdamped. We observe, though, that the average rate of 
energy loss is quite similar, as indicated in figure |21 where we plot the energy that exists in the form of particles (the 
last term in (|2.24|) ') as a fraction of the conserved total energy. In fact, this figure shows that the overall conversion 
of inflaton energy into particles is somewhat more efficient in the non-equilibrium evolution than is suggested by the 
adiabatic approximation. 
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FIG. 3: Energy of particles as a fraction of the conserved total energy. 



Of course, the question arises, what accounts for the large qualitative difference between the non-equilibrium 
evolution and that obtained from the adiabatic approximation? A simple answer is that, as shown in section ITT1 
no adiabatic expansion of the non-equilibrium evolution equations actually exists, and our numerical results merely 
confirm that the adiabatic approximation is wrong. This may perhaps be the appropriate answer, but it needs some 
further discussion, because we also showed in section [H] that the adiabatic expansion does exist if we take the extra 
step of replacing self-energies with their equilibrium values, and one might think that this step is relatively innocuous 
if the state is not too far from equilibrium. 

Consider first the initial state, which was chosen to be a state of exact thermal equilibrium. According to the 
adiabatic approximation, the criterion for neglecting cf> is very well satisfied (equation 1)3.80 ). an d the motion should 
be strongly overdamped. However, we see from figure |21 that \(f>\ as calculated from the non-equilibrium evolution 
equations is very large. The reason for this is not hard to discover. In a state of exact thermal equilibrium, the 
term g(j)(x 2 (x, i))trunc in the equation of motion 1|2.I9[1 is precisely the quantity AV'((j>) that appears in the adiabatic 
equation 12.2511 . Thus, while the original equation of motion asserts that U(<fi) = V'(cf>) + AV'(4>) is equal to — <j), 
the adiabatic approximation asserts that U{4>) ~ —i]((j>)(j), the value of <f> being very small. Clearly, the adiabatic 
approximation gives a wrong answer for a state of exact thermal equilibrium. In fact, the informal arguments given 
m 0- El an d reviewed in |l5j| show that the frictional effects arise precisely from small departures from equilibrium. 

Figure^lillustrates the departure from thermal equilibrium at times near t 20, corresponding to the first negative 
peak in figure|2 where about one half-cycle of oscillation has been completed. The quantity plotted is k 2 uiknk, whose 
integral gives the total energy in particles. The solid curve uses the distribution rik(t) generated by the non-equilibrium 
evolution, while the broken curve uses the equilibrium distribution ro^ q with the temperature /3 _1 adjusted to give 
the same total energy. (This is quite independent of the equilibrium distribution n% dmh generated by the adiabatic 
evolution.) Clearly, these two distributions are different, but the difference does not greatly affect the time evolution. 
For example, we show in figure [S] the quasiparticle width Tfc, which appears in the evolution equations (|2.5|) and 
(|2.tJ|) . calculated from rik(t) and n£. q , and the difference is seen to be marginal. We can use the effective value of (3 
calculated as above, together with the non-equilibrium values of and Tk in the expressions l|2.f 8fl and (|2.22|) to find 
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FIG. 4: Weighted particle distribution k 2 ujkrik at time t = 19.9 (solid curve) and the same quantity calculated from the 
equilibrium distribution for the same total particle energy (broken curve). 

an effective, non-equilibrium value for the ratio [x 2 /rj 2 which, according to (|3.5[) . gives a criterion for overdamping at 
times near t = 20. The result is that \fJ?/r] 2 \ « 9.5 x 10~ 5 . Thus, the criterion for overdamping (as obtained from 
the adiabatic approximation) is well satisfied by the non-equilibrium state at t « 20, just as it was by the initial 
equilibrium state. 

The main point of the preceding discussion is this. In order to obtain the adiabatic approximation, we had to 
replace non-equilibrium self-energies Sfe[n] with the values Sfe[n cq ] associated with a state of thermal equilibrium. In 
the initial state at t = these two self energies are identical, and in the state illustrated in figure the difference 
is quite small. Clearly, this difference does not account for the drastic failure of the adiabatic approximation in the 
situation we have examined. The difference between the two evolutions in figure [21 (and, for that matter, those shown 
in figure^) indicates that the expectation value (x 2 (x, i))trunc is not well represented as the sum of the two terms in 
(|2.21|) . regardless of whether the coefficients AV'((j>) and r]((f>) are calculated using the non-equilibrium state or an 
effective equilibrium state. Thus, although the rate of thermalization may well be an important factor, its importance 
does not lie in the relatively small effect that inefficient thermalization would have on the values of these coefficients. 

Having made this point, we defer further discussion to section llVl but we observe here that thermalization in the 
situation we have studied is rather less efficient than one might infer from the value 0.97 for the ratio in l|3.6[l . In the 
initial state, the relaxation time r k = T^ 1 has values w 400 for all the modes used in the calculation, and this is 
longer than the whole time for which the evolution was studied. From the point of view of the adiabatic approximation, 
this means that although thermalization might have been expected to be moderately efficient by comparison with the 
initially slow evolution, it is much less efficient by comparison with the more rapid evolution at later times. In the 
non-equilibrium evolution, there are quite large departures from equilibrium when cj) changes rapidly. For example, 
we show in figureElthe particle distribution (again weighted by k 2 uik) at the time t w 10.2 where 4> first becomes very 
small. This figure illustrates the fact that the particles created by the source (|2.12|l have a non-thermal spectrum 
weighted towards high energies: it falls off roughly as uj^ 5 compared with the exponential fall-off of the thermal 
distribution. It is on account of this difference that the integral equation (|2.15l) fails to respect the energy sum rule, 
and thus to be self-consistent. Comparing figures 0] and it appears that, although the created particles arc not 
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FIG. 5: The quasiparticle width F^ calculated from the particle distributions of figure [I] 



thcrmalized instantaneously, a significant amount of thermalization has taken place during the quarter-cycle between 
t = 10.2 and t = 19.9. This appearance is somewhat illusory, however, as we explain in section llVl 



C. Numerical details 



The calculations reported above were performed using modest computing resources (a PC), and a moderate numer- 
ical accuracy seems to be adequate for the issues we seek to address. The most demanding aspect of the calculations 
is the solution of the evolution equations l|2.5|) and i|2.6[l . To deal with these, we computed the evolution of N modes 
(that is, N values of the momentum k) for the distribution functions n/- and v^, with an upper cutoff /c max = 200. Of 
the various integrals that are needed, the energy integral l|2.24|l attaches the largest weight to high-momentum modes 
and the chosen cutoff is large enough to make this integral reasonably accurate, as illustrated by the integrands in 
figures 0] and Since the calculations are quite lengthy (the data needed to produce figure were obtained from 
3 weeks continuous running) we did not attempt to make the entire process converge by increasing fc max and N 
systematically. However, the results presented here, using N = 400, which adequately resolves the structure of the 
distributions generated, differ negligibly from those obtained using N — 200, and seem to be insensitive to the cutoff 

^max ■ 

It is apparent from l|2.6|) that the functions t%(t) oscillate with frequencies roughly equal to u>k- The basic time 
stepper used was a semi-implicit method, which is stable for these oscillations, with a fixed time step St about one- 
tenth of the period of the fastest oscillation. However, computation of the scattering integrals (|2.7I) and l|2.8|) is very 
time-consuming (even after analytical reduction to triple integrals), so ak and T^, which change quite slowly relative 
to the fast oscillations, were updated at longer time intervals At, using a predictor-corrector method that typically 
converged after two or three iterations. This larger time interval At was assigned values between 205t and 150<5i, 
according to the rate of change of Tfe (see figure 0). The values of rik needed in these integrals were obtained by 
interpolating the N tabulated values with a cubic spline. 

The quality of the data obtained by these methods is illustrated in figured which shows the variation over a typical 
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FIG. 6: Weighted particle distribution at t = 10.2, compared with the equilibrium distribution for the same energy. 



cycle of 4> and the k = modes of Tk and rik- The sharp peaks in Tk appear in all modes, and are largely due to the 
factor 1/oJk in l|2.17|) . The basic operations of evaluating integrals, and so on, were carried out with a tolerance of 
0.5% or better. The energy l|2.24|l . which formally is exactly conserved, exhibits small variations, but is constant to 
within about 1%, which we estimate as the overall accuracy of the calculation. 



IV. SUMMARY AND DISCUSSION 



According to quantum field theory, the equation of motion for the expectation value of a scalar field is a complicated, 
non-local integro-differential equation. In applications to inflationary cosmology, it is typically assumed that this 
equation of motion can be replaced by a local equation similar to l|l.l|) . If it is further assumed that the state of the 
quantum system is fairly close to thermal equilibrium, then it seems plausible that one should be able to apply the 
methods of linear response theory, which yield quantitative estimates of the friction coefficient /?((/>) and the effective 
potential V e s((j>). 

Here and in 0, EU we have tried to address the question whether the local equation of motion obtained in this 
way adequately approximates the actual non-equilibrium evolution, at least when this evolution is sufficiently slow. 
The indications are that it does not, but the treatment we have given cannot be considered as definitive, and we wish 
to set out clearly the issues that are involved. 

The first issue is that we have treated the non-equilibrium evolution on the basis of an approximate set of local 
equations l|2.2[) . (|2.4|) . (|2.5|l and 1)2.6(1 . As described in detail in they arise from a local Ansatz for quasiparticle 
self-energies, and it is not easy to assess how good an approximation this is. For the purposes of estimating the non- 
equilibrium behaviour numerically, other approximation schemes have recently become available (see, for example, 
[la. ITa Il8l fl9 | ) which are in some respects superior to this one. However for the purpose of deriving a local effective 
equation of motion, approximations more or less equivalent to ours would seem to be inevitable. To be precise, the 
question we address here is whether the fully local equation of motion H2.25f) as derived from linear response theory 
yields a good approximation to the non-equilibrium behaviour as described by our evolution equations. If it does not 
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FIG. 7: Variation of cj>, nk = o an d r^ = o over a typical cycle of oscillation. Data points are separated by the time interval At 
explained in the text. 



(as we seem to find), then it is also unlikely to provide a good approximation to the true non-equilibrium dynamics. 

The phrase "fully local" above means the following. The evolution equations (|2.2I) . (|2.5|) and (|2.6|) are local, in 
the sense that they involve 4>(t), rik(t) and Vk{t) together with their time derivatives, all evaluated at the same time 
t: they describe a Markovian process. In principle, one can envisage solving (|2.5|l and (|2.6|l for n,k(t) and Vkit) as 
functionals of <fi(t). These solutions are non-local: they depend on the entire history of <f>{t). So, substituting these 
solutions into the remaining equation l|2.2|l . we have a non-local equation of motion for <j). This non-local equation 
might be approximated by the "fully local" equation l|2.25|l if the quantity gifiix 2 ) nas a time-derivative expansion 
whose first two terms are AV'((j>) and rj{(j))(j). Although we cannot find the complete functionals rik{t), and Vk(t), 
we can attempt to find their time-derivative expansions. As explained in section [HJ these expansions do not exist. 
However, if we make the further approximation of calculating self-energies from an equilibrium particle distribution 
rather than the non-equilibrium distribution, we can obtain a time derivative expansion, which essentially reproduces 
the results of linear response theory. This is what we have referred to as the adiabatic approximation. In section ITTT1 
we compared numerically the non-equilibrium evolution with that resulting from the adiabatic approximation, and 
found very large differences. We also saw that, for the model studied here, thcrmalization is not very fast compared 
with the motion of (f>. Intuitively, it seems that this poor thcrmalization might account for the failure of the adiabatic 
approximation. 

The second issue we want to discuss is ; therefore, that of thermalization. As indicated in section lTlI Bl the relaxation 
time tx, = 1/Tfc is roughly equal to 4>/<fi in the initial state from which figure |21 was generated. On the basis of the 
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adiabatic approximation, which indicates strong overdamping, one would perhaps expect the non-equilibrium evolution 
to exhibit a period of overdamped motion, albeit with significant corrections due to a modest thermalization rate. 
What we actually found was oscillatory motion with a period much shorter than r k . Over one quarter-period of this 
motion, the particle distribution changes from that shown in figure 0] to that shown in figure [fj] From these figures, 
it would seem that significant thermalization has taken place over a period of time much shorter than r k . However, 
figure suggests a different explanation. From that figure, we see that particles with a non-thermal distribution 
are radiated over a short period of time as 4> 2 decreases towards 0, and that these particles, rather than being 
thermalized, are quickly reabsorbed as 4> 2 subsequently increases. So thermalization is indeed slow on the time scale 
of the non-equilibrium motion. 

We cannot altogether rule out the possibility that it is this slow thermalization that accounts for the failure of 
the adiabatic approximation in the situation we have studied. However, we do not think that this is the correct 
explanation. One reason is that, as we have just remarked, consideration of the adiabatic approximation alone would 
suggest that thermalization, while modest, is not so slow as to produce the gross discrepancy apparent in figure [21 A 
more compelling reason is that, as discussed in section llll Bl the non-equilibrium motion is not well represented by the 
adiabatic equation (|2.25J) even at those times when the state is close to thermal equilibrium. This strongly indicates 
that the step of replacing, say, T k [n] with Ffc[n eq ] is wrong in principle and not just because of the quantitative 
discrepancy which, as seen in figure [SJ may be quite small. In more detail, suppose that the distributions n k (t) and 
v k (t) are indeed close to thermal equilibrium. Then it should be a good approximation to linearize the evolution 
equations i|2.5fl and l|2.6l) . using n k (t) = n c k q (t) + Sn k (t) and ffe(i) = f° q (£) + 5v k {t). In particular, the linearized 
version of 12.511 is 

poo ■ 

d t Sn k + / dk'K(k, k')8n k , = d t n q - —Re + 8v k \. (4.1) 
Jo u k 

Unlike (|2.15|) . this equation is self-consistent; the energy sum rule merely gives an expression for the rate of change of 
the quantity J dk k 2 aj k 5n k . (This equation is not self-contained, however. To solve it, we would need a prescription, 
such as the conservation of energy, for determining the time-dependent temperature that specifies n cq . As a matter of 
fact, the scattering integral S k vanishes when n cq has an arbitrary chemical potential. In the adiabatic calculations of 
section ITTT1 we did not allow for the possibility of a non-zero chemical potential because there is no conserved particle 
number that could be used to determine its value.) The important point here is that the linearized equation i|4.1|) 
is a differential equation. Its solution is again a non-local functional of cj>, which will not reproduce the adiabatic 
equation l|2.25f) . regardless of how small Sn k might be. To do that, we need not just a linear approximation, but a 
time-derivative expansion which, as we have emphasized, is not self-consistent. 

In summary, our numerical study indicates that that the non-equilibrium evolution of (j) is not well described by a 
local equation of motion with the friction term estimated by linear response theory. It is possible that an adiabatic 
approximation would work better for a model in which thermalization is more efficient, and we plan to explore this in 
future work. However, for the reasons just discussed, we do not think that poor thermalization is the principal cause 
of the large discrepancies we observed. 
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